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ABSTRACT 

We have obtained the first continuous disk averaged spectrum of Mars from 450 to 1550 Ghz using the Herschel-SPTRE Fourier 
Transform Spectrometer. The spectrum was obtained at a constant resolution of 1.4 GHz across the whole band. The flux from the 
planet is such that the instrument was operated in "bright source" mode to prevent saturation of the detectors. This was the first 
successful use of this mode and in this work we describe the method used for observing Mars together with a detailed discussion of 
the data reduction techniques required to calibrate the spectrum. We discuss the calibration accuracy obtained and describe the 
first comparison with surface and atmospheric models. In addition to a direct photometric measurement of the planet the spectrum 
contains the characteristic transitions of 12 CO from J 5-4 to J 13-12 as well as numerous H2O transitions. Together these allow 
the comparison to global atmospheric models allowing the mean mixing ratios of water and 12 CO to be investigated. We find that 
it is possible to match the observed depth of the absorption features in the spectrum with a fixed water mixing ratio of 1 x 10~ 4 
and a 12 CO mixing ratio of 9 x 10 -4 . 
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1. Introduction 

Water vapour plays a key role in the Martian atmospheric 
chemistry and physics. As source of hydrogen radicals it 
has an important impact on the ozone (Lefevre et al., 
2004) concentration, the stability of the Martian atmo- 
sphere (Parkinson and Hunten, I1972p and plays an im- 
portant role in the comparative planetology context (e.g. 
Hartogh et al. 120041 120 10|) . The water cycle and the asso- 
ciated variable hygropause (Clancy 119961 Encrenaz [2001) 
with a strong asymmetry of its distribution at perihelion 
and aphelion is one of the principal problems of the present 
Martian climate. If confirmed, the asymmetry is expected 
to drive a strong unbalance in the meridional transport. 
A better observational characterization will be a strong 
constraint for general circulation models of the Martian 
atmosphere and particular its meridional transport (e.g. 
Medvedev and Hartogh. I2007p . Since its first detection in 
the near infrared (Spinrad et al, 1963) water vapour in 
the Martian atmosphere has been extensively studied by 
ground-based and space borne observations (e.g. Sprague, 
etaLfTOOOl Burgd orf, et al. 120001 Gurw ell, et al. 120001 Smith 
120021 Biver, et a l 120051 Fouchet. 120071 Encrenaz et al. 120081 
and Smith et al, 2009 ) . Compared to water, there is only a 
sparse database for carbon monoxide at Mars. Space borne 



* Herschel is an ESA space observatory with science instru- 
ments provided by European-led Principal Investigator consor- 
tia and with participation from NASA. 



( Encrenaz . 120061 and Smith. 12009P and ground-based (e.g. 
Lellouch et al, 119911 Krasnopolsky, I2007P show an aver- 
age abundance of 700 ppm. One of the major goals of the 
Herschel is to trace the origin, evolution and chemistry of 
water in the Solar system (Hartogh et al 2009) using all 
three Herschel instruments, including SPIRE, to take spec- 
tra of a variety of Solar system objects including Mars. 
Mars was thought to be too bright to be observed using the 
SPIRE Fourier-transform spectrometer (FTS) as it was an- 
ticipated that its bolometric detectors would saturate and 
provide no useful information. However, as described below, 
we have devised a new mode of operation for the instrument 
that prevents saturation and which has allowed us to use 
Herschel-SPIKEi to take the first continuous spectrum in 
the submillimetre of the Martian atmosphere and provide 
the first complete set of water vapour and carbon monoxide 
spectra in the 450 to 1550 GHz range. Given the novelty 
of this observing mode and the difficulties encountered in 
reducing and calibrating the data, this paper describes in 
some detail the observing mode, data reduction and calibra- 
tion before briefly describing an initial Martian atmospheric 
model used to validate our calibration approach. 



2. Observations and observing Mode 

Mars was observed using the central pixels on the SPIRE 
spectrometer on Herschel operational day 176, 2009 
November 6. See Griffin et al 120101 for an explanation of 
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Fig. 1. The observed submillimetre spectrum of Mars reduced and calibrated as described in the text compared to a 
model spectrum generated using the RADTRANS code (Irwin |2009[) . The positions of the detected H2O and CO lines 
are indicated. 



the SPIRE instrument and its observing modes. The ob- 
servation period commenced at 19:36:16 and finished at 
20:20:14 UT: a total of 1332 seconds were spent on the 
spectrum. During the observation the target was tracked 
by the Herschel satellite in non-sidereal pointing mode. 
The target was at RA 8h 51m 37.3s and Dec 19°15'40.5" 
at the start of the observation. The spectrum was taken 
as part of the functional and calibration set up observa- 
tions used to establish the optimum operating conditions 
for the instrument. At this time in the mission the notion- 
ally 2% emissivity spectrometer calibration source (SCAL- 
2) was set to 15 K to assess its efficacy in "nulling" the flux 
from the telescope and thus providing an accurate refer- 
ence signal against which the source spectrum is measured. 
Subsequent observations showed that use of the SCAL is 
not necessary and, in fact, adds photon noise to the spec- 
trum without adding to the photometric accuracy. The sub- 
millimetre flux density of Mars is extremely high in the 
context of observing with the SPIRE instrument which was 
designed to cope with flux densities up to a few hundred Jy 
rather than the approximately 10000 Jy from Mars. Under 
nominal detector bias and sensitivity conditions the power 
from Mars would saturate the detectors and conditioning 
electronics and no useful data would result. However, the 
method of biasing and reading out the bolometers is to use 
sinusoidal AC bias voltage with a 160 Hz frequency. The sig- 
nals from the bolometers are then converted to a DC volt- 
age through a lock-in amplifier that derives its frequency 
reference from the bias supply electronics. The long cable 
(5 m) between the warm electronics and the bolometers in 
combination with the resistance of the bolometers induces 
a phase error between the AC bolometer voltage at the in- 
put to the lock-in amplifier and the frequency reference. In 



normal operations this is removed by calibrating this phase 
error and maximising the bolometer signal whilst observing 
"dark sky" . However, when observing very bright objects 
it is desirable to reduce both the inherent sensitivity of the 
bolometers and the overall electronics gain. No gain con- 
trol was built into the SPIRE conditioning electronics so 
instead we have found that we can vary the gain by delib- 
erately de-phasing the reference and bias waveforms. The 
variation in gain between a square wave demodulation and 
a sine wave is proportional to the cosine of the phase error 
between them. In this way then we have a sensitive elec- 
tronics gain control that can be used to attempt to avoid 
saturation when observing "bright" sources. 

3. Data processing and calibration 

For the Mars observation we set the out of phase condition 
between the reference and bolometer bias to give a notional 
gain of about 0.4 for SSW and about 0.28 for SLW com- 
pared to the gain when observing dark sky. In practice, 
however, the electronics gain is only one part of the overall 
change in response of the detector system when observ- 
ing a very bright source such as Mars. We find that even 
after applying the known variation in electronics gain the 
derived flux density of Mars is still much lower than that 
expected from any model continuum. During all observa- 
tions with the SPIRE spectrometer the PCAL source (see 
Griffin et al 120101 and Swinyard et al l2010[) is operated to 
give a fixed and repeatable power level onto the detectors. 
Inspection of the variation in recorded bolometer temper- 
ature between the PCAL flash during a reference observa- 
tion of Uranus and that during the Mars observation after 
correction for the electronics gain variation shows that the 
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Table 1. Principal line identifications in the Mars FTS 
spectrum. The line to continuum ratios for H2O lines 11 
and 12 are shown in Fig. 2 and for CO line 15 in Fig. 3. 
*This line is blended with 11 



Line Number 


Species 


Transition 


Frequency (GHz) 


1 


Ortho H 2 


lio-loi 


556.94 


2 


CO 


J=6-5 


691.491 


3 


Para H 2 


2n-2o2 


752.03 


4 


CO 


J=7-6 


806.65 


5 


CO 


J=8-7 


921.8 


6 


Para H 2 


2o2-ln 


987.93 


7 


CO 


J=9-8 


1036.91 


8 


Ortho H 2 


3l2"3o3 


1097.37 


9 


Para H 2 


I11-O00 


1113.34 


10 


CO 


J=10-9 


1152.01* 


11 


Ortho H 2 


3l2"221 


1153.13 


12 


Ortho H 2 


321-3l2 


1162.91 


13 


Para H 2 


4 2 2-4l3 


1207.64 


14 


Para H 2 


220-2n 


1228.79 


15 


CO 


J=ll-10 


1267.01 


16 


H 2 


625-532 


1322.06 


17 


CO 


J=12-ll 


1382.0 


18 


H 2 


523"5l4 


1410.6 



nal result is shown in Fig. 1 together with the model spec- 
trum. We can see that there is now good agreement in the 
flux density across the SPIRE range but some distortion 
of the continuum spectrum in the SSW range as well as a 
small scaling difference between the overall continuum and 
the measured spectrum amounting to between 5 and 10% 
across waveband. Inspection of the voltage timelines shows 
that, even with the reduced electronics gain there is still 
a degree of saturation in SSW when observing Mars and 
this causes the low order part of the transformation to be 
distorted. The level of agreement is satisfactory and, im- 
portantly for the overall interpretation of the spectrum, is 
demonstrated to be only subject to scaling and not offset 
errors. The absolute calibration accuracy using this method 
can be estimated from the uncertainties in the ratio of the 
PCAL flashes and to first order are around 20% for the 
SLW detector - similar to those quoted for the nominal op- 
erational mode of the SPIRE FTS (Swinyard et al l2TJTO)) . 
However for the SSW we do not have an independent mea- 
surement of the uncertainty and, given that there are also 
distortions evident in the continuum spectrum, we can as- 
sign an accuracy of no better than 30% at present. New 
observational and calibration techniques for bright sources 
are under development and we expect this to improve. 



overall response of the system was a factor of 1.5 lower dur- 
ing the Mars observation for the central detector in SLW. 
The signal to noise and signal stability for the PCAL data 
in SSW were not sufficient to allow any estimate of the 
variation in response and a different approach was adopted 
here (see below). In order to use this information for the 
Mars data we have processed the interferograms through 
an alternative processing scheme compared to the standard 
one described in Fulton et al (2008). Here we use an obser- 
vation of Uranus taken on Herschel operational day 209 
(2009 December 8) and a model by Moreno (fT998l f20T0|) to 
provide the calibration reference. We convert the averaged 
voltage interferograms from the detectors into net bolome- 
ter temperature for the dark sky, Mars and Uranus ob- 
servations separately. The Uranus and associated dark sky 
were transformed into spectral space, the dark sky spectra 
subtracted from the Uranus observation and the net spec- 
trum divided by the Uranus continuum model to provide 
the calibration relative spectral response function (RSRF). 
For the Mars observation the SLW net temperature inter- 
ferograms were first multiplied by the estimated response 
variation deduced from the PCAL flashes (xl.5) and then 
transformed into the spectral domain . The dark sky ob- 
servation for Mars was taken with nominal phase settings 
and no response factor was applied before transformation. 
The dark sky was subtracted from the Mars spectrum and 
the net spectrum divided by the RSRF. The same proce- 
dure was followed for the SSW data but initially with no 
gain factor applied to the Mars interferogram. Inspection 
of the result in comparison to a model (described in Sect. 
4) showed a good agreement between the SLW spectrum 
and the model but a significant variation with the SSW 
spectrum. We then applied a variable scaling factor to the 
SSW Mars interferogram until the recovered flux density 
in the overlapping spectral region between SSW and SLW 
came into agreement. The scaling factor finally applied was 
1.35, in reasonable accord with that used for SLW. The fi- 



4. Atmospheric model of Mars 

The observation was modelled using the RADTRANS 
line-by-line radiative transfer code (Irwin [2009) and the 
HITRAN spectral database. We prescribed an input atmo- 
spheric temperature profile and a surface brightness tem- 
perature by using outputs from the European Martian 
Climate Database (EMCD, v4.1). The EMCD is a database 
describing the climate and environment of the Martian at- 
mosphere, derived from martian general circulation models 
(Forget 119991 Lewis et al. ll999)l . Essentially, for the appear- 
ance and season (L s , sol) of Mars relevant to our observa- 
tion, we partitioned the visible martian disk in a fine, reg- 
ular grid (typically 100 points across a martian diameter). 
The EMCD database was then used to provide the local 
surface pressure and the vertical profiles of temperature in 
the atmosphere and in the sub-surface. The surface pres- 
sure and the temperature-pressure-altitude profiles were av- 
eraged over the disk, providing the input thermal profile 
needed for the radiative transfer calculation. For calculat- 
ing the surface brightness temperature, radiative transfer 
was performed within the surface, using the online tool of 
Lellouch et al. (|2008|) : a surface brightness temperature of 
220 K was derived. Finally, the derived spectrum was con- 
volved with a sine function of FWHM 0.048cm" 1 and aver- 
aged over the disk. For the CO lines we tried mixing ratios 
between 4.5 and 18xl0" 4 and for the H2O lines we tried 
mixing ratios between 0.5 xlO -4 and 2xl0 -4 . The disk- 
averaged flux density was then finally calculated given the 
Mars solid angle of 1.25xl0~ 10 for the observation date of 
2009 November 6. 

The full spectrum model generated with fixed mixing 
ratios of 9xl0" 4 for CO and lxlO" 4 for H 2 is plotted 
together with the calibrated observed spectrum in Fig. 1 
and the identified lines are given in Table 1. Note that the 
CO 4-3 and 4-5 lines were not clearly detected in the spec- 
trum as these regions suffered from instrumcntally induced 
spectral fringes which masked the absorption features. The 
level of agreement between the model continuum level and 
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Fig. 2. Close up of the measured and model spectrum 
around the Ortho 3i2-22i and 32i-3i2 water lines showing 
the effect of varying the mixing ratio (see text) 
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Fig. 3. Close up of the measured and model spectrum 
around the CO J=ll-10 line showing the effect of varying 
the mixing ratio (see text) 



the observed spectrum is good except in the SSW range 
where some distortion is found for reasons described above. 
A more detailed view of the H2O line to continuum ra- 
tios is shown in Fig. 2. We have taken three model cases 
here to show the variation of line depth and width with 
constant mixing ratios of 2xl0 -4 (purple), lxlO -4 (green) 
and 0.5xl0 -4 (red). It is evident that the green line shows 
the best fit to the line depths. However, the line wings in the 
model appear broader than in the observation and further 
more detailed modelling in conjunction with higher spectral 
resolution H erschel-iUFl observations will be required to 
fully match the spectrum. In Fig. 3 we show a close up 
of the spectrum around the CO J=ll-10 line showing the 
effect of CO mixing ratios of 18xlCP 4 (purple), 9xl0~ 4 
(green) and 4.5xl0~ 4 (red). Over all CO lines 9xl0 -4 was 
found to be the best fit. 



before launch and it represents a significant enhancement 
to the observation capabilities of Herschel. We have taken 
the first continuous submillimetre spectrum of the planet 
Mars and shown the ability to measure all water and car- 
bon monoxide lines between 600 and 1550 GHz with good 
photometric accuracy in a single short observation. Our 
preliminary analysis of the spectrum shows that a model 
with abundances of 900 ppm for CO and 100 ppm for H2O 
best fits the observed line depths. Both these values are in 
line with previously established figures. A more sophiscated 
modeling approach is clearly needed to make the best use 
of these and future data sets that will be obtained during 
the Herschel mission and this will be the subject of future 
work. 
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5. Conclusions 

This paper has demonstrated the utility of the SPIRE FTS 
in making efficient broad band spectral observations of very 
bright astronomical sources and the ability to provide ac- 
curately calibrated spectra when SPIRE is operated in its 
"bright source" mode. This was not expected to be the case 



STFC Rutherford Appleton Laboratory, Harwell Innovation 

Campus, Didcot, OX11 0QX, U.K. 

Max Planck Institute for Solar System Research, Katlenburg- 

Lindau D-37191, Germany 

Bluesky Spectroscopy, Lethbridge, Canada 

School of Physics and Astronomy, Cardiff University, Cardiff, 

UK 



B. M. Swinyard et al.: The Herschel-SPTKE submillimetre spectrum of Mars 

5 LESIA, Observatoire de Paris, France 

6 Physikalisches Institut, University of Bern, Switzerland 

7 Space Research Centre, Polish Academy of Sciences, Warsaw, 
Poland 

LERMA, Observatoire de Paris, France 
9 European Space Astronomy Centre, European Space Agency, 
Madrid, Spain 

10 Instituto de Astrofisica de Andalucia, Granada, Spain 

11 California Institute of Technology, Pasadena, USA 

12 Instituut for Sterrcnkunde, Katholieke Universiteit Leuven, 
Belgium 

13 Argelander Institute for Astronomy, University of Bonn, 
Germany 

14 Astronomy Department, University of Michigan, USA 

15 Universit de Bordeaux, Laboratoire d'Astrophysique de 
Bordeaux, France 

16 Laboratory of Molecular Astrophysics, CAB-CSIC. INTA, 
Madrid, Spain 

17 Joint Astronomy Center, Hilo, USA 

18 Joint ALMA Observatory, Chile 

19 Institute d'Astrophysique et de Geophysique, Universit de 
Lige, Belgium 

20 Jet Propulsion Laboratory, California Institute of 
Technology, Pasadena, USA 

21 Department of Physics and Astronomy, University of 
Lethbridge, Canada 

22 1st Physics Institute, University of Cologne, Germany 

23 SRON, Netherlands Institute for Space Research, 
Netherlands 

24 Physikalisches Institut, University of Bern, Switzerland 



